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FIELD OF THE INVENTION 

[01] Our invention relates generally to cellular wireless networks. More particularly, our 
invention relates to methods for estimating non-uniform spatial offered loads generated 
throughout a service territory associated with the network. 

DESCRIPTION OF THE BACKGROUND 
[02] Today, cellular wireless carriers primarily use second-generation cellular systems to 
provide cellular communications. The most widely used systems are the global system 
for mobile communications (GSM), which uses a hierarchy of time division multiplexing 
access (TDMA) frames, and code division multiple access (CDMA), which is a spread 
spectrum-based technique. In addition, significant work is underway on the development 
of third-generation wireless communications systems. 

[03] Regardless of the technology used, cellular carriers provide wireless services by 
partitioning their service territories into cells. Ideally, a cell is in the shape of a hexagon; 
however, in practice, topological limitations and other considerations lead to significant 
variations in a cell's topology. Each cell may be further partitioned into a few (e.g., 3 or 
6) cell sectors. Every cell sector has a Base Transceiver Station (BTS) with multiple 
transceivers that transmit and receive signals at multiple frequencies. Thus, the BTS 
serves as a mobile station's access point into the communications network. 

[04] Designing and operating a wireless cellular network is quite complex. One approach a 
carrier typically uses to mange its network is to partition its service territory into small 
"bins" (also called pixels) and to evaluate the service quality at each bin. For example, a 
50-kilometer by 50-kilometer service territory can be partitioned into more than 60,000 
bins of 200 meters by 200 meters. Using propagation models and field measurements, a 
carrier then derives the mean and standard deviation of the signal strength received at 
every bin from every BTS. The matrix of the expected signal strengths among all bin- 
BTS combinations is known as the Received Signal Strength Indicator (RSSI) matrix. A 
carrier uses this matrix to determine if a mobile station in a particular bin is being 
adequately served. In particular, multiple BTS's may adequately serve a mobile station 
in a particular bin. Typically (e.g., in GSM), the BTS that transmits the strongest signal to 
the mobile location serves the mobile station. Alternatively (e.g., in CDMA), multiple 
BTS's that transmit the strongest signals to the mobile location may serve the mobile 
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station simultaneously. A mobile station in a specific bin is served well if the ratio of the 
carrier signal strength to that of the sum of the interfering signal strengths is above a 
certain threshold. 

[05] A carrier also continuously collects performance data at every BTS, including data 
regarding the carried load, blocked calls, dropped calls, quality of connections, etc. From 
this information, the carrier attempts to infer the service provided to different bins in the 
territory. 

[06] However, carriers today do not know the offered load (represented, for example, in units 
of bits per second, number of time slots per second, minutes of call per minute, etc.) 
generated at individual bins at different hours of the day. Knowing each bin's offered 
load would better facilitate the carrier's planning and operational processes. In particular, 
knowing the offered load at each bin would allow the carrier to better manage network 
resources and to provide better service to those bins with higher loads by allowing more 
intelligent assignment of frequencies among the BTS's and by facilitating load balancing 
across the BTS's through the fine-tuning of different control parameters. 

[07] More specifically, a fundamental problem with many of today's cellular network 
technologies (e.g., in GSM) is the assignment of frequencies among the BTS's. The total 
number of available frequencies is limited so that each frequency must be assigned to 
multiple BTS's. However, the assignment of a given frequency across multiple BTS's 
must be done so that the mobile stations experience negligible interference. For example, 
Karen I. Aardal, Stan P. M. Van Hoesel, Arie M. C A. Koster, Carlo Mannino, and 
Antonio Sassano present in the paper "Models and Solution Techniques for Frequency 
Assignment Problems," that appeared as a report of Zentrum fur Informationstechnik 
Berlin (ZIB), ZIB Report 01-40, December 2001 a survey on frequency assignment 
models and algorithms, which paper is hereby incorporated by reference. As the authors 
describe, the models/algorithms use an aggregation of the bins to cell sectors and use as 
input a representative expected signal strength received at every BTS from every other 
BTS. However, frequency assignment models and algorithms would be improved by 
assigning weights to the bins wherein these weights are set equal to the bin estimated 
offered loads. In addition, Jean-Marie Bourjolly, Leslie Dejoie, Ke Ding, Oumar 
Dioume, and Michel Lominy emphasize in their paper "Canadian Telecom Makes the 
Right Call, Frequency Allocation in Cellular Phone Networks: an OR Success Story," 
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OR/MS Today, 29, 40-44, April 2002 that the resulting frequency plan should be 
evaluated at the bin level rather than at the BTS level with the objective of minimizing, 
for example, the number of bins with unacceptable reception quality. Knowing the bins 
that generate higher offered loads would allow a carrier to spend more resources and 
thereby provide better service to bins that generate more demand. 

[08] BTS load management is also an important issue in order to avoid uneven congestion and 
blocking of call attempts. Carriers can, for example, adjust the transmission power from 
specific BTS's and thus change the area served by each of the BTS's. By reducing a 
certain BTS's transmission power, some of the bins that this BTS previously served will 
now receive a stronger signal from other BTS's, which will now serve those bins. Thus, 
the fine-tuning of the power parameters at a BTS is used to balance loads among the 
BTS's. Any load-balancing scheme that takes into consideration the loads generated at 
each bin would more effectively fine-tune the power parameters. 

SUMMARY OF OUR INVENTION 
[09] Accordingly, it is desirable to provide methods that overcome the shortcomings of the 
prior art and estimate non-uniform spatial offered loads over a wireless service territory 
of interest in a cellular carrier's network. Our method uses as input a partition of the 
service territory into a plurality of logical bins, the mean and standard deviation of the 
signal strength received at every bin from every BTS in the service territory, an offered 
load estimate for each Base Station Transceiver (BTS) that serves the service territory, 
and an approximation of relative demand for wireless services across the bins in the 
service territory. As output, our method provides estimates of the offered loads for each 
of the bins in the service territory of interest. 

[10] More specifically, the mean and standard deviation of the signal strength received at 
every bin from each of the BTS' are first determined using measurements and known 
propagation models. Treating the received signal strength at a given bin from a given 
BTS as a random variable, probability density functions and cumulative distribution 
functions are defined for each bin across each BTS using the determined means and 
standard deviations. Based on the defined functions, probabilities are derived for each 
bin across each BTS, indicating the probability a given BTS will serve a given bin. Each 
of these probabilities is then computed. 
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[11] Next, an estimate of the offered load for each BTS is determined using measurements of 
carried load and lost calls at the BTS. Using these BTS offered loads and the computed 
probabilities, a resource constraint is expressed for each BTS in terms of the bin offered 
loads that our method determines. 

[12] The last input to our method is an approximation of the relative demand for wireless 
services across the bins in the service territory. These relative demands are obtained from 
demographic data, such as information on residence density, business density, and traffic 
patterns, and are obtained independently of BTS load measurements. These relative 
demands are converted to demand targets such that the sum of the demand targets across 
all the bins equals the sum of the BTS estimated offered loads across all the BTS's. Using 
these demand targets, a performance function is then defined for each bin in terms of the 
bin's offered load that is to be determined. 

[13] Once having the BTS resource constraints and the performance functions, an equitable 
resource allocation model is generated, wherein the model is formulated as a 
lexicographic minimax optimization problem with decision variables that represent the 
estimated non-uniform spatial offered loads for each bin. Specifically, a lexographic 
minimax objective function is defined in terms of lexicographic minimization of the 
vector of the performance functions sorted in a non-increasing order. This objective 
function is then coupled with the resource constraints to define the model. Using known 
algorithms, the model is finally solved thereby finding the set of estimated bin offered 
loads that both results in the lexographic smallest vector of performance functions sorted 
in a non-increasing order and that satisfies the resource constraints. 

[14] Significantly, a carrier can use the estimated bin offered loads as bin weights for 
performing frequency assignment across all BTS's. Also, a carrier can use the estimated 
bin offered loads to identify the bins with higher offered loads. Knowledge of the bins 
with higher offered loads allows a carrier to improve the service provided to these bins. 
Furthermore, a carrier can use this knowledge of the bin offered loads to balance the 
loads imposed on the BTS's in the service territory in a more effective way, for example, 
by a more effective adjustment of the signals power parameters. 

BRIEF DESCRIPTION OF THE DRAWINGS 
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[15] Figure 1 depicts the method steps of one illustrative embodiment of our invention for 
estimating the non-uniform spatial offered loads over the bins of a cellular wireless 
service territory. 

DETAILED DESCRIPTION OF OUR INVENTION 
[16] Figure 1 is a flow chart depicting the method steps of our invention for estimating the 
non-uniform spatial offered loads over a service territory of interest in a carrier's cellular 
wireless network over a specific time period. Beginning with step 102, the carrier's 
service territory, or a portion of this territory of interest, is divided into logical bins. For 
description purposes, these bins can be numbered 1 through N (i.e., N is the number of 
bins in the service territory of interest). The exact number of bins and the geographic size 
of each bin are not specific to our invention. 

[17] A first input to our method is the mean and standard deviation of the signal strength 
received at every bin from every BTS in the service territory of interest. For description 
purposes, the BTS's in the service territory of interest can be numbered 1 through M (i.e., 
M is the number of BTS's in the service territory of interest). Accordingly, in step 104, 
physical test measurements are made at a sample of the N bins wherein each 
measurement is the received signal strength from a sample of the M BTS's. These sample 
measurements are correspondingly used to infer the received signal strength among all 
bin-BTS combinations. Using these measurements and known propagation models (see, 
for example, Chapter 10 of William Stallings, "Wireless Communications and 
Networking", Prentice Hall, 2002 for propagation models), in step 106 the mean //*, and 
standard deviation a X j of the received signal strength at each bin j for all j = 1 to N across 
each BTS i for all for i=l to M is then determined. 

[18] The received signal strength at any given bin j (for j = 1 to AO from any given BTS i (for 
i = 1 to M) is a random variable and can be represented by a probability density function 
fij{s). As such, proceeding to step 108, using the mean fig and standard deviation oft of 
the received signal strength, a probability density function fi/s) is defined for each bin j 
for all j = 1 to N across each BTS i for all i = 1 to M, where frfs) is the density function 
for the signal strength s that bin j receives from BTS /. Note that in wireless propagation 
models, signal strengths are typically expressed in dBm (decibel-miliWatt), which is a 
logarithmic scale. As such, the signal strengths are typically modeled as random 
variables with log-normal density functions. Accordingly, when the signal strength is 
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expressed in dBm, each f^s) is typically a normal distribution with mean fa and standard 
deviation (Ty. Nonetheless, it should be noted that our method is independent of the 
actual density function assumed for the signal strength. Having f^s), the cumulative 
distribution function Fy(5) can also be defined for all j = 1 to N across each BTS i for all 
for i=l to M, where Fy(5) is the probability that the signal strength received at bin i from 
BTS j does not exceed S. 

[19] Proceeding to step 110,/^) and F 0 (5) are used to compute p tj for each bin j for all j = 1 
to N across each BTS i for all i=l to M, where p u represents the probability that BTS i 
will serve bin j or, in other words, represents the probability that a mobile station in bin j 
is served by BTS L In particular, the derivation of p,y is based on assumptions underlying 
the wireless technology covering the service territory of interest. Accordingly, as these 
assumptions change, /?,;, should be modified. For example, in multiple technologies, such 
as GSM, it is reasonable to assume that the BTS that provides the strongest received 
signal to a given bin serves any given mobile station within that bin and that each of the 
signal strengths received from the M BTS's within the bin are independent random 
variables. From these assumptions, p t j can be expressed as shown in equation (1): 

Pij = ifij(s)[ n F k j(s)]ds (1) 
s Ic * i 

for /= 1, 2,. . ., M and j = 1, 2,. ... N. In practice, we set probabilities p$ to zero if the mean 
signal strength received at bin j from BTS i is, for example, 18 decibels below the largest 
mean signal received at bin j (which implies that the signal is, on average, only about 
1.5% of the strongest signal). Once defining the probabilities py as expressed in equation 
(1), for example, the probabilities are computed using numerical integration methods. 
Finally, the positive probabilities for each bin j are normalized to a sum of one. 

[20] As a second example, in cellular wireless technologies such as CDMA systems, the up to 
C (where C > 1) BTS's that provide the strongest received signals to a given mobile 
station in a given bin can simultaneously serve that mobile station. Here, the probability 
giji must be determined, wherein is the probability that bin j receives its /-th strongest 
signal from BTS i for all i= 1 to M,j = 1 to N, and I = 1 to C Accordingly, let B(j) be the 
subset of BTS's that are the candidate BTS's for serving the mobile stations in bin j (e.g., 
B(j) may include all BTS's that produce mean signals to bin j that are within 18 decibels 
of the strongest mean signal.). In addition, let B denote a subset of {B(j) - i) with I - 1 
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elements. Then equation (2a) represents the probability for each bin j across each BTS i 
that the strongest signal bin j receives comes from BTS i: 

Si/l = J///W[ n F k j(s)]ds (2a) 
s Vk,k±i 

for i= 1, 2,..., M and j = 1,2,..., N. In addition, equation (2b) represents the probability 
for each bin ; across each BTS i for each of the remaining C - 1 strongest signals that this 
signal at bin j comes from BTS i: 

g iH = 1 Ifij(s) n [1-FhWI n F k j(s)ds (2b) 

B:BQB(j)s keB kzB(j)-B-i 

fori= 1,2,..., MJ= 1,2,..., N, and/ = 2, 3, ...,C. Once g m and g# for / = 2, 3, ...,Care 
derived, the probabilities py can be expressed as shown in equation (3): 

C 

Pij= 1 Sijl (3) 
J 1 = 1 J 

for /= 1, 2,. . ., M, and j= l,2,...,N. Once defining the probabilities g I>7 as expressed, for 
example, in equations (2a)-(2b), the probabilities are computed using numerical 
integration methods. In practice, we set probabilities py to zero if the mean signal strength 
received at bin j from BTS i is, for example, 18 decibels below the largest mean signal 
received at bin / Finally, the positive probabilities py for each bin j are normalized to a 
sum of one. 

[21] A second input to our method is an offered load estimate for each BTS in the service 
territory of interest, which loads are determined in step 112. In particular, the offered 
load estimate L, is estimated for each BTS i for all i =1 to M using carried load 
measurements and information regarding lost calls at each BTS i over the time period of 
interest (note that the carried load is defined as the offered load multiplied by the 
proportion of calls that are not lost). 

[22] Using these offered load estimates L, and the probabilities p (j determined in step 110, a 
resource constraint is expressed in step 1 14 for each BTS i for all / =1 to M. In particular, 
let Xj represent the estimated offered load at bin j for all j = 1 to N (i.e., xj for all j=ltoN 
are the non-uniform spatial offered loads our method ultimately estimates). The offered 

8 



App 1563 



load for a given BTS i can be computed using the estimated bin offered loads xj and the 
probabilities p tj as shown in expression (4): 



for i = 1, 2, M. Importantly, the computed offered load at a given BTS / as shown in 
expression (4) cannot exceed the estimated offered load L { for BTS i as determined in 
step 112. Accordingly, inequality (5) represents a resource constraint expression for any 
given BTS /for? =1,2, ...,M: 



In other words, inequalities (5) are constraints on the non-uniform spatial estimated 
offered loads xj at bin j for all j = 1 to N that our method ultimately estimates. 
Specifically, the set of constraints for i =1, 2, . . M ensures that all estimated bin offered 
loads are such that the total computed offered load at any specific BTS i (as a function of 
the estimated bin offered loads) does not exceed that BTS's offered load estimate L,. 

[23] A third input to our method is an approximation of relative demands for wireless services 
across the bins in the service territory of interest, which demands are determined in step 
116. In particular, a relative demand estimate is determined for each bin j for all j = 1 to 
N. As an example, demographic data corresponding to the time period of interest can be 
used to determine the relative demand estimates for wireless services. Such data includes 
information on residence density, business density, and traffic patterns. Internal 
marketing organizations and wireless services industry consultants can provide these 
approximations. Note that the relative demands are independent of the load 
measurements at the BTS's. 

[24] Proceeding to step 118, the relative demand estimate for each bin j for all j = 1 to N is 
converted to a demand target e y In particular, each relative demand estimate is 
proportionally adjusted so that the sum of the demands over all N bins is equal to the sum 
of the measured offered loads L, over all M BTS's as shown in equation (6). 



N 



(4) 




J ' = l 



IsPijXj^Li 



(5) 



N M 



Y,ej = El, 



(6) 
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[25] In step 120, the demand targets <? 7 are used to express a performance function h/xj) for 
each bin ; for all j = 1 to N. As shown by equation (7), each performance function h/xj) 
for a given bin j measures the weighted normalized deviation of the bin's estimated 
offered load Xj from the bin's demand target ej. 

hj(xj) = Wj— (7) 

for all; = 1, 2 N and where the weights wj > 0. Note that hpcj) is strictly decreasing 

with xj. In most cases, each of the weights vv, are set to one, which implies that the 
performance functions represent normalized deviations from the demand targets. When 
Wj = e jf for all j, the performance functions represent deviations from demand targets. Our 
method can handle any selected weights. 

[26] Once having the resource constraint expression for each of the M BTS's as determined in 
step 114 and as shown in inequality (5) and the performance function h/xj) as determined 
in step 120 and as shown in equation (7), an equitable resource allocation model is 
generated, wherein the model is formulated as a lexicographic minimax optimization 
problem with decision variables that represent the estimated non-uniform spatial offered 
loads xj for all j = 1 to N for the bins in the service territory of interest. The optimal 
solution to the model provides the offered load estimate Xj for each of the N bins within 
the service area. 

[27] Specifically, proceeding to step 122, the performance functions h/xj) for each of the N 
bins are used to express a lexicographic minimax objective function, which, when 
coupled with the resource constraints (5), form an equitable resource allocation model. 
Specifically, this model is defined as an allocation of "available resources" (i.e., offered 
load estimates at the BTS's) such that no performance function value h n (x n ) for a given 
bin n can be decreased further without either violating at least one of the resource 
constraints (5) or without increasing the performance function value hk(x k ) of another bin 
k wherein h^x k ) is at least as large as h n (x n ) for all n and k = 1 to N. Since the 
performance functions are strictly decreasing, an equivalent way of defining the equitable 
resource allocation model is that no estimated offered load x n for a given bin n can be 
increased without either violating at least one of the resource constraints (5) or without 
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decreasing another estimated offered load x k at a given bin k that has a performance 
function value hk(x k ) that is at least as large as h n (x n ). 

[28] Accordingly, let x be a vector of the N estimated offered loads Xj and let h N (x) be a vector 
of the N performance functions hj{xj) for j = 1 to N, where these performance functions 
are sorted in non-increasing order as shown in formulation (8a) and (8b). 

h N (x) = [h h (x h ), h h (x h ),..., h jN (x jN )] (8a) 

where, 

h h {x h )>h h {x h )>...>h jN (x jN ) (8b) 

The lexicographic minimax objective function of the equitable resource allocation model 
consists of finding the lexicographic smallest vector h N (x). Letting x l be the vector x > 0 
of estimated non-uniform spatial offered loads across the N bins that results in the 
lexicographic smallest vector h N (x) among all vectors x > 0, the lexicographic minimax 
objective function can be expressed as shown in equation (9). 

^( JC / ) = lexmin[^W] (9) 

X 

among all vectors x > 0. 

[29] Proceeding to step 124, the equitable resource allocation model is expressed by coupling 
the resource constraints (5) with the lexicographic minimax objective function of 
equation (9). This model is formulated by (10a) - (10c). 

^U / ) = lexmin[^(jc)] (10a) 

X 

i Pij Xj<Li for/=l,2,...,M (10b) 

xj>0 fori = 1,2, ...,Mand;= 1,2, (10c) 

Objective function (10a) minimizes lexicographically the vector h N (x), where hj{xj) is 
given by equation (7) and h N (x) is defined by (8a)-(8b). Constraints (10b) ensure that the 
sum of the estimated offered loads at all bins, multiplied by the appropriate probabilities, 
does not exceed the offered load estimated at each of the M BTS's. Constraints (10c) 
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ensure that all estimated loads are nonnegative. The optimal solution to the model as 
formulated by (10a)-(10c) is such that no performance function value h/xj) can be 
feasibly decreased without increasing the value of another performance function value 
that is at least as large. 

[30] The elements of the vector x l that provides the optimal solution to the model 

formulated by (10a)-(10c) are the non-uniform spatial offered load estimates jc y for each 
bin j for j = 1 to N within the service territory of interest over the time period of interest. 
The vector x l is an equitable allocation of the offered load estimates at the BTS's among 
the N bins, since no estimated offered load x n for a given bin n can be feasibly increased 
without decreasing another estimated offered load jc* at a given bin k that has a 
performance function value hk(x k ) that is at least as large as h n (x n ). 

[31] Proceeding to step 126, the optimal solution x l to the equitable resource allocation model 
formulated by (10a)-(10c) is determined, thereby providing the offered load estimates xj 
for each of the N bins within the service territory of interest. In particular, Hanan Luss in, 
"On Equitable Resource Allocation Problems: A Lexicographic Minimax Approach," 
Operations Research, 47, pages 361-378, 1999 (hereinafter "Luss"), presents in Section 2 
several variations of algorithms for solving the model as formulated by (10a)-(10c) using 
today's desktop computers, which paper is hereby incorporated by reference. In addition, 
Hanan Luss and Andrew Vakhutinsky in, "A Resource Allocation Approach for the 
Generation of Service-Dependent Demand Matrices for Communications Networks," 
Telecommunication Systems, 17, pages 411-433, 2001, present an adaptation of an 
algorithm to solve a model for determining service-dependent point-to-point demands in 
wire-line networks and further present computational results for large networks. 
Nevertheless, it should be noted that our invention is independent of the exact algorithm 
used to solve the equitable resource allocation model as formulated by (10a)-(10c). 

[32] In general, an algorithm for finding the optimal solution to the equitable resource 
allocation model, as formulated by (10a)-(10c), repeatedly solves minimax models. Upon 
each minimax solution, some variables are fixed at their lexicographic minimax values 
and excluded from the model, and the amount of available resources is updated. 
Eventually, all variables are fixed at their lexicographic minimax values, thus 
representing the equitable resource allocation solution. 
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[33] More specifically, let / denote a set of BTS's and / denote a set of bins. Initially, / = { 1, 
2,..., M} and J = {1, 2,..., N). We assume that the input parameters to our method are 
cleaned-up so that U > 0 for all i e /, so that for each i e /, pg > 0 for some j e J t and so 
that for each j e J, p {j > 0 for some i e I. The initial minimax model to be solved is 
formulated by (llaMllc): 

V* = rmn[maxMx,)] ( lla ) 

x j 

IPijXjZLi, is I (lib) 



Xj 



>0, jeJ (11c) 



It is important to note that the model as formulated by (lla)-(llc) may have multiple 
optimal solutions. Accordingly, as part of an algorithm for solving the model as 
formulated by (llaMllc), we need to compute the minimal minimax solution vector**. 
Note that jc* is the "minimal minimax solution" to the model formulated by (1 la)-(l lc) if 
jc* < jc\ where x* is any other minimax solution to the model formulated by (1 la)-(l lc). 

[34] What follows is an outline of one variation of the above-mentioned algorithms for 
solving the equitable resource allocation model as formulated by (10a)-(10c). As 
indicated, these algorithms include the step of computing the "minimal minimax 
solution" vector x* for the model as formulated by (lla)-(llc). Accordingly, an outline 
of one variation of an algorithm for computing this vector x* is first presented. Note that 
this particular algorithm for computing the "minimal minimax solution" vector x* is 
described in "Luss" and in Hanan Luss and Donald R. Smith, "Resource Allocation 
Among competing Activities: A Lexicographic Minimax approach," Operations 
Research letters, 5, pages 227-231, 1986, which paper is also hereby incorporated by 
reference. Again note that our invention is independent of the exact algorithm used to 
solve the equitable resource allocation model as formulated by (10a)-(10c) and of the 
exact algorithm used to compute the "minimal minimax solution" for the model 
formulated by (1 la)-(l lc). 
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Exemplary Algorithm for Computing the "Minimal Minimax Solution" Vector x* 
for the Model as Formulated by (lla)-(llc). 

1 . Initialize temporary sets JTEMP = J and ITEMP = /. 

2. Compute equation (12). 



X Pijej-Li 

Vj= JeJTEMP $ . e [TEMp 

£ Pijej/wj 

je JTEMP 



(12) 



3. Find the solution to the relaxed minimax problem (i.e., without constraints (1 lc)) 
as shown in equation (13): 

V R =maxlVA (13) 

ie ITEMP 

4. If wj > V* for all j e JTEMP, then proceed to Step (6) below; otherwise, exclude 
from JTEMP any j e JTEMP for which vv, < V*. 

5. Update expressions V, as given by equation (12), for all i e ITEMP, by 
subtracting from the summations those terms that correspond to the /s that were 
deleted from JTEMP in Step (4). Delete any constraint is ITEMP for which all ; 
with pij > 0 were deleted from JTEMP. Return to Step (3). 

6. Record the minimal minimax solution as shown by equations (14a)-(14b) and 
stop. 

V* = V R (14a) 
^ = max[/ l 7H^*),0] = max[e y (l-V*/M; ; ),0], je J (14b) 

Note that if the w/s are all the same, then, in Step (4), wj > V* for all je JTEMP in the 
first iteration. In this case, Step (4) implies that the algorithm finds the minimal minimax 
solution in a single iteration. 

Exemplary Algorithm for Solving the Equitable Resource Allocation Model as 
Formulated by (10a)-(10c). 

1. Formulate the minimax model as represented by equations (lla)-(llc) with sets / 
and J and loads L v 
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2. Compute the minimal minimax solution to the model represented by equations 
(lla)-(llc) using the above described six step "minimal minimax solution" 
algorithm (or any other algorithm that provides the minimal minimax solution). 
The minimal minimax solution is given by equations (14a)-(14b). Identify the set 
RO ={i: is ITEMP and V, = V*}. These are resource constraints that are fully 
used. 

3. Identify set JFIX = {j: p xj > 0 for some / s RO}. Fix x l j = x) for all j e JFIX. 

4. Delete from J all j e JFIX. If J is now empty, stop; the equitable resource 
allocation solution has been obtained. If J is not empty, continue with Step (5). 

5. Delete i from / if all j with /?y > 0 were deleted from /. Update the loads 
Li <— Li - £ Pu xj for all i 6 /. Return to Step (2). 

j&JFIX 

[36] As indicated, the bin offered load estimates can be used as weights for the bins wherein 
these weights would be used to improve frequency allocation models and algorithms. 
Also, once having the offered load estimate for each bin, a carrier can identify the bins 
with higher offered loads. Knowledge of the bins with higher offered loads allows a 
carrier to improve the service provided to these bins by modifying the frequency 
assignment plan across the BTS's. Furthermore, a carrier can use this knowledge of the 
bin offered loads to balance the loads imposed on the BTS's in the service territory in a 
more effective way. 

[37] As an example of our method, consider a cellular network service territory of interest 
divided into N = 6 bins with M = 3 BTS's. Assume the normalized probabilities, Py, 
determined from steps 104-1 10 for each of the six bins across each of the three BTS's is 
as shown in Table 1: 



1 BTSl 


BTS2 


BTS3 


Binl 


0.5 


0.5 


0 


Bin 2 


0 


0.5 


0.5 


Bin 3 


0.5 


0 


0.5 


Bin 4 


0.2 


0 


0.8 


Bin 5 


0.2 


0.8 


0 


Bin 6 


0 


0.2 


0.8 



Table 1 



15 



App 1563 



[38] Assume further that the determined offered loads L, as determined in step 1 12 for each of 
the three BTS's are 100, 200, and 300 respectively. Accordingly, using inequalities (5) 
and as determined in step 114, the resource constraint for each of the three BTS's is as 
follows: 

BTS 1: 0.5jc, + 0.5jc 3 + 0.2jc 4 + 0.2x 5 < 100 

BTS 2: 0.5JC! + 0.5jc 2 + 0.8jc 5 + 0.2jc 6 < 200 

BTS 3: 0.5jc 2 + 0.5jc 3 + 0.8jc 4 + 0.&t 6 < 300 

where xj > 0 for j = 1 to 6. 

[39] Assume also that the demand targets ej for each of the six bins and as determined in steps 
116 and 118 are 150, 100, 50, 150, 100, and 50 respectively. Accordingly, using weights 
Wj = 1 for all j = 1 to 6, the performance functions for each bin and as determined in step 
120 are given by equation (7). Using these performance functions to express a 
lexicographic minimax objective function (i.e., step 122) and coupling this objective 
function with the above resource constraints, an equitable allocation resource model for 
the bin offered loads is expressed (i.e., step 124). Solving this model in step 126, the 
non-uniform spatial offered load for each bin is as shown in Table 2. 



Binl 


*!=ioo 


Bin 2 


4 = 161.111 


Bin 3 


4 = 33.333 


Bin 4 


4=ioo 


Bin5 


4 = 66.667 


Bin 6 


4 = 80.556 



Table 2 



Note that this solution satisfies the resource constraints for BTS's 1 and 2 as equality and 
has an excess of 58.333 of the resource associated with BTS 3. 

[40] Comparing each bin's offered load to its demand target, it is seen that bins 1, 3, 4, and 5 
have estimated offered loads below the demand target, while bins 2 and 6 have estimated 
offered loads above the demand targets. 
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[41] As another example, assume the determined offered loads L, as determined in step 112 
for each of the three BTS's are each 200 at a different time-of-day (the remaining 
assumptions regarding p xi and e 7 are the same as above). Here, the solution to the 
equitable allocation resource model for the non-uniform spatial offered load for each bin 
is as shown in Table 3. 



Binl 


jc{ = 144.135 


Bin 2 


4-85.106 


Bin 3 


4 = 42.553 


Bin 4 


4 = 127.660 


Bin 5 


4 = 96.090 


Bin 6 


4 = 42.553 



Table 3 



The solution satisfies the resource constraints for BTS's 2 and 3 as equality and has an 
excess of 61.906 of the resource associated with BTS 1. 

[42] Again, comparing each bin's offered load to its demand target it is seen that the estimated 
offered loads for all six bins are now below their demand targets. Importantly, the 
solution of our model as formulated by (10a)-(10c) in this second example provides 
estimated offered loads at the bins that are different from the estimate offered loads of the 
first example since our model considers the changes in the offered load estimates at each 
of the BTS's, which are derived from measurements collected at the BTS's. The demand 
targets, on the other hand, are only rough estimates that may depend, for example, on 
demographic data and on the sum of offered load estimates over all the BTS's. However, 
the demand targets are derived without using the information of offered load estimates at 
specific BTS's. 

[43] The above-described embodiments of our invention are intended to be illustrative only. 
Numerous other embodiments may be devised by those skilled in the art without 
departing from the spirit and scope of our invention. 
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